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Helioseismology and solar neutrinos: an update 
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We review recent advances concerning helioseismology, solar models and solar neutrinos. Particularly we address 
the following points: i) helioseismic tests of recent SSMs; ii) predictions of the Beryllium neutrino flux based on 
helioseismology; iii) helioseismic tests regarding the screening of nuclear reactions in the Sun. 



1. Introduction 

Helioseismology has added important data on 
the solar structure which provide severe con- 
straint and tests of SSM calculations. For in- 
stance, helioseismology accurately determines the 
depth of the convective zone Rb and the photo- 
spheric helium abundance Y v h- With these ad- 
ditional constraints there are essentially no free 
parameters for SSM builders. 

In this paper we review recent advances con- 
cerning helioseismology, solar models and so- 
lar neutrinos. Particularly we shall address the 
following points: i) helioseismic tests of recent 
SSMs; ii) prediction of the Beryllium neutrino 
flux based on helioseismology; iii) helioseismic 
tests concerning the screening of nuclear reactions 
in the Sun. 

2. A summary of helioseismic determina- 
tions of solar properties 

While we refer to e.g. Q| for a review of 
the method and to Q for the data, we recall 
that by measurements of thousands of solar fre- 
quencies (p-modes) with a typical accuracy of 
Au/u ~ 10 -4 , one derives: 

• properties of the present convective enve- 
lope, such as depth and helium abundance: 




R b = 0.711(1 ±O.4%)i? 
Y ph = 0.249(1 ±4%) 



(1) 
(2) 



The quoted errors, mostly resulting from 
systematic uncertainties in the inversion 



Figure 1. The isothermal sound speed profile, u = 
P/p, as derived from helioseismic observations, 
from §. 



technique, have been estimated conserva- 
tively by adding linearly all known individ- 
ual uncertainties, see 0. If uncertainties 
are added in quadrature, the global error 
is about one third of that indicated in eqs. 
(§]J|), see again [0. This latter procedure 
was also used by Bahcall et al. (|] with sim- 
ilar results. This yields the so called "lcr" 
errors. We shall refer to the conservative 
estimate as the "3cr" determination. We re- 
mark however that this terminology is part 
of a slang, and it does not correspond to well 
defined confidence levels, as one has to com- 
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bine several essentially systematic errors. 

• sound speed profile. By inversion of hclio- 
seismic data one can determine the sound 
speed in the solar interior. This analysis can 
be performed in terms of either the isother- 
mal sound speed squared, u = P/p, or in 
terms of the adiabatic sound speed squared 
c 2 = dP/dp\ad — jP/p, as the coefficient 
7 = d\ogP/dlogp\adiab is extremely well 
determined by means of the equation of 
state of the stellar plasma. 

In fig. [l] we show the helioseismic value of u 
as a function of the radial coordinate R/ Rq. 
The typical "3er" errors are of order ±0.4% 
in the intermediate solar region, R/ Rq ~ 
0.4, and increase up to ±2% near the solar 
center. 

3. Helioseismic tests of recent Standard 
Solar Models 

Fig. H compares the results of six different ob- 
servational determinations of the sound speed 
with the results of the best solar model of ref. 
Q, hereafter BP2000. This figure suggests sev- 
eral comments: 

i) Different measurements yield quite consistent 
value of the sound speed, to the level 0.1%; 

ii) The solar model BP2000 is in agreeement 
with helioseismic data to better than 0.5% at 
any depth in the sun. We remark that also the 
properties of the convective envelope predicted by 
BP2000 {R b /R Q = 0.714, Y ph = 0.244) are in 
agreement with helioseismic determinations, see 
eqs. §|). 

iii) On the other hand, the predicted sound 
speed differs from the helioseismic determination 
at the level of 0.3-0.4% just below the convective 
envelope. 

Concerning this last point, we remark that the 
difference is however at the level of the "3cr" un- 
certainty of the helioseismic determination, as it 
is shown in Fig. [| Nevertheless it can be taken as 
an indication of some imperfection of the SSM. A 
marginally better agreement can be obtained in 
solar models including mixing induced by rota- 
tion, see e.g. B and in models including macro- 
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Figure 2. Fractional difference between the sound 
speeds calculated for the BP2000 model and the 
sound speeds in six helioseismological experi- 
ments ||], from Q). 



scopic transport term, see e.g. ||. 

In summary all this means that SSM predic- 
tions are accurate to the level of one per cent or 
better, although there are indications of some de- 
ficiencies at the level of per mille. 

Concernig neutrino physics, the properties of 
the central solar region, where nuclear reactions 
are efficient, are relevant. Specifically, Boron and 
Beryllium neutrinos are produced very near to 
the solar center with maximal production rates 
respectively at Rb — 0.04i? Q and Rb & = O.O6i?0. 
Since the p-modes which are observed do not 
propagate (actually are exponentially damped) 
so deeply in the sun the question often arises if 
present helioseismic data can determine the sound 
speed in region of Beryllium and Boron produc- 
tion. From an extensive analysis of the inversion 
method and of data available at that time we al- 
ready concluded in jy that u(R ~ 0) is deter- 
mined with a "Ict - " accuracy of 1%. This point 
was further elucidated in || where a simplified 
analysis was presented in order to produce con- 
vincing evidence that helioseismology fixes the 
sound speed near the solar center with such an 
accuracy. 
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Figure 3. Comparison between BP2000 Q and 
MDI helioseismic data for u = P/p. The "lcr" 
and "3fi" helioseismic uncertainties jlj corre- 
spond to the dark and light areas respectively. 



4. Helioseismology and the beryllium neu- 
trino flux 

As well known, the production of neutrinos 
from 7 Be + e~ —^ 7 Li + ^ e is an important item in 
the context of the so called "solar neutrino puz- 
zle" for several reasons: 

i) The result of Gallium experiments would be 
(partially) consistent with the hypothesis of stan- 
dard neutrinos if Beryllium neutrino (vse) pro- 
duction rate, L(vbc), is suppressed by an order 
of magnitude with respect to the prediction of 
the Standard Solar Model (SSM), see e.g. @,|]. 

ii) If one accepts neutrino oscillations as the solu- 
tion of the solar neutrino puzzle, the determina- 
tion of the neutrino mass matrix depends however 
on the predicted value of L{vB e )- 

m) Direct experiments aiming to the determina- 
tion of the VBe signal are in preparation and of 
course the interpretation of their result will rely 
on L{v Be ) |Jl0|. 

The SSM prediction for Beryllium neutrinos 



i, 

L{v B e) SSM = 1-3 • 10 37 (1 ± 9%) s- 1 (la) (3) 

is very robust, much more than that of Boron 
neutrinos. However any additional information 
which docs not rely on SSM are clearly welcome. 

In ref. ]ll|] a lower limit on the Beryllium neu- 
trino flux on earth was found, Q(Be) m i n — 1 • 10 9 
cm" 2 s _1 , in the absence of oscillations, by using 
helioseismic data, the B-neutrino flux measured 
by Superkamiokande and the hydrogen abun- 
dance at the solar center X c predicted by Stan- 
dard Solar Model (SSM) calculations. We remark 
that this abundance is the only result of SSMs 
needed for getting &(Be) m i n . Lower bounds for 
the Gallium signal, G m j„ = (91 ±3) SNU, and for 
the Chlorine signal, C mm = (3.24 ±0.14) SNU, in 
the absence of oscillations have also been derived. 
They are about 3<r above the corresponding ex- 
perimental values, G exp = (75 ±5) SNU 
and C exp = (2.56 ± 0.22) SNU |j. We remark 
that predictions for X c are very stable among dif- 
ferent (standard and non standard) solar models, 
see ftl|| . In fact X c is essentially an indicator of 
how much hydrogen has been burnt so far. The 
stability of X c corresponds to the fact that any 
solar model has to account for the same present 
and time integrated solar luminosity. 

In ref. Jl5|] a step forward was made; a determi- 
nation of L(i>Be) was directly obtained by means 
of helioseismology without using additional as- 
sumptions. The basic idea is the following. In the 
SSM the pp-II termination (which is the Beryl- 
lium producting branch) accounts for an appre- 
ciable fraction of the 4 He produced near the solar 
center. If Beryllium production is suppressed - 
now and in the past - less He would have been 
produced near the center. As a consequence, the 
molecular weight should decrease and the sound 
speed should increase in this region. In other 
words, we know that SSM calculations are in good 
agreement with helioseismology and we can ex- 
pect that this agreement is spoiled if i>Be produc- 
tion is substantially altered. 

In order to test this idea, solar models with 
artifically changed vse production were con- 
structed. An efficient way for producing arbi- 
trary variations of L{vse) is to vary the zero en- 



Figure 4. The photospheric helium abundance 
Y P h and the depth of the convective envelope 
Rb in solar models with the indicated values of 
s = S34/ S^f M . The error bars correspond to the 
"lcr" helioseismic uncertainties, from ref.flt]. 



Figure 5. Fractional difference with respect to 
the SSM prediction, (model-SSM)/SSM, of the 
isothermal squared sound speed, u = P/p, in 
solar models with the indicated values of s = 
Sa4,/S^ M . The dotted area corresponds to the 
"lcr" helioseismic uncertainty on u, from ref.H. 



ergy astrophysical S-factor (S34) for the reaction 
3 He+ 4 Hc Be+7, since, as well known (see e.g. 
]l^||jl7j ) , the production rate Liyse) is directly 
proportional to S34: 

L(Be)/L(Be) SSM ~ S 3i /Sl£ M . (4) 

We remind that S34 is measured with an accu- 
racy of about ten per cent, (S§f M = 0.54 ± 0.05 
KeVb |18|[), It has been varied, however, well 
beyond its experimental uncertainty in order to 
simulate several effects which have been claimed 
to suppress VBe production, e.g. hypothetical 
plasma effects which could alter nuclear reaction 
rates. The resulting values for the quantities 
which can be tested by means of helioseismology 
are shown in Fig. ^ and Fig. for several assumed 
values of the parameter s = S34/ S^f AI . 

The photospheric helium abundance Y p h is 
weekly sensitive to the value of S34 whereas the 
depth of the convective envelope Rb is altered by 
more than lcr if S34 is reduced below one half of 
the SSM value, i.e.: 

L{vbc) = 1-3 ■ 10 37 (1 ± 50%) s- 1 at lcr . (5) 



sound speed is altered, particularly near the solar 
center. In fact, stringent constraints arise from 
the sound speed profile, particularly near R ~ 
O.2i? . The requirement that the sound speed is 
not changed by more than la yields: 

L{vBe) = 1-3 • 10 37 (1 ± 0.25) s- 1 at lcr . (6) 

In conclusion, helioseismology directly confirms 
the production rate of Beryllium neutrinos as pre- 
dicted by SSMs to within ±25% (lcr error). This 
constraint is somehow weaker than that estimated 
from uncertainties of the SSM, see Eq.||, however 
it relies on direct observational data. 

5. Helioseismology and screening of nu- 
clear reactions in the sun 

The study of screened nuclear reaction rates 
was started with the pioneering work of Salpeter 
jfj|, who discussed both the extreme cases of 
"weak" and "strong" screening, providing suit- 
able expressions for the screening factors 



As previously mentioned, one expects that the 



fij ij , plasma / {cu} ijfiare 



(7) 
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The solar core is not far from the weak screening 
case, however it does not satisfy the usual condi- 
tions under which the weak screening approxima- 
tion holds. This is the reason why the possibilities 
of large deviations from weak screening have been 
investigated by several authors, expecially as an 
attemp to avoid or mitigate the "solar neutrino 
problem" , see e.g. 
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and references therein. 
We remind that solar models are built by us- 
ing stellar evolutionary codes which include spe- 
cific expressions for the nuclear reaction rates. If 
one uses different formulas for the screening fac- 
tors fij one obtains different solar models. On 
the other hand, helioseismology provides precise 
information on the sound speed profile and on 
the properties of the convective envelope which 
have to be reproduced by the correct solar model. 
Therefore helioseismology can potentially provide 
a test of screening models. 

In ref. [glj , solar models using different screen- 
ing assumptions were constructed and compared 
with helioseismic data. Specifically, four different 
model were considered: 

i) The weak screening approximation (WES). The 
screening factors fij are given by: 



]nf^ ES = Z i Z j e 2 /{a D kT) 



(8) 



where Zi,Zj are the charges of the interacting 
nuclei, T is the temperature and an is the Debye 
radius. As clear from equation above, the screen- 
ing factors are always larger than unity, i.e. the 
plasma provides enhancement of the thermonu- 
clear reaction rates. 

ii) The Mitler result (22) (MIT), obtained with 
an analytical method which goes beyond the lin- 
earized approach and which correctly reproduces 
both the limits of weak and strong screening. 

iii) Neglect completely any screening effect 
(NOS), i.e. nuclear reactions occur with rates 
(ct>)bare- This case is considered in connection 
with the suggestions that screening can be much 
smaller than Salpeter's estimate, see e.g. p3| . 

iv) The Tsytovich model (TSY) f|§|||], which 
provides a decrease of all the thermonuclear reac- 
tion rates with respect to the case of bare nuclei. 

In Table | the screening factors at the solar 
center for the various models are shown. One 
sees that the weak screening approximation al- 



Table 1 

Screening factors in solar center, for weak screen- 
ing (WES) @, Mitler model (MIT) ||, no 
screening (NOS) and Tsytovitch model (TSY) 

M. 





WES 


MIT 


NOS 


TSY 


p + p 


1.049 


1.045 


1 


0.949 


3 He + 3 He 


1.213 


1.176 


1 


0.814 


3 He + 4 He 


1.213 


1.176 


1 


0.810 


7 Be + p 


1.213 


1.171 


1 


0.542 



ways yields the largest enhancement factors, as 
physically clear due to the fact that electrons and 
ions are assumed to be free and capable of follow- 
ing the reacting nuclei. By definition there is no 
enhancement in the NOS model, whereas in TSY 
model there is a decrease of the reaction rate, as 
already remarked. 

Recent Standard Solar Models calculated by 
using the weak screening prescription are in 
agreement with helioseismic constraints on the 
properties of the convective envelope and on the 
sound speed profile, see fig. || This shows that the 
weak screening model is in agreeement with data 
and deviations from WES cannot be too large. 
From the comparison of different models, see Ta- 
ble | and fig. | one obtains the following results: 

i) The difference between the Tsytovitch model 
(TSY) and the weak screening model (WES) ex- 
ceeds the "conservative" uncertainty on u in a 
significant portion of the solar profile. We re- 
mark that also the depth of the convective enve- 
lope is significantly altered. In other words the 
anti-screening predictions of ref. @,^5| can be 
excluded by means of helioseismology. 

ii) Also the difference between the no-screening 
model (NOS) and WES is significant for both u 
and Rb in comparison with helioseismic uncer- 
tainty. In other words the existence of a screening 
effect can be proved by means of helioseismology. 

iii) The Mitler model of screening (MIT) can- 
not be distinguished from the weak screening 
model within the present accuracy of helioseis- 
mology. 
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Figure 6. Comparison of different screening mod- 
els with the WES model for u = P/p. The 
"lex" and "3c" helioseismic uncertainties Q cor- 
respond to the dark and light areas. 



Table 2 

Comparison among solar models with different 
screening factors. We show the fractional differ- 
ences, (model -WES )/WES, for the photospheric 
helium abundance (Y p h), depth of the convective 
envelope (Rb), central temperature (T c ), isother- 
mal sound speed squared at the solar center (u c ), 
neutrino fluxes (<&i) and predicted signals for 
Chlorine (CI) and the Gallium (Ga) experiments. 
All variations are in per cent. 





MIT 


NOS 


TSY 


Yph 


-0.076 


-0.86 


-1.4 


R b 


+ 0.037 


+0.34 


+0.59 


T c 


+ 0.45 


+0.54 


+1.4 


u c 


+ 0.10 


+1.0 


+1.4 


%p 


+0.033 


+0.45 


-0.35 




-0.19 


-2.4 


-5.9 




-2.7 


-12. 


+11 


CI 


-2.5 


-11. 


+9.7 


Ga 


-0.76 


-2.9 


+2.3 



6. Concluding remarks 

We summarize here that main points of our dis- 
cussion: 

i) The comparison of recent SSMs with helioseis- 
mic measurements shows that SSMs prediction 
are accurate to the level of one per cent or better, 
although there are indications of some deficiencies 
at the level of per millc. 

ii) Helioseismology directly confirms the produc- 
tion rate of Beryllium neutrinos as predicted by 
SSMs to within ±25% (la error). 

iii) Models for screening of nuclear reactions in 
the Sun can be tested by means of helioseismol- 
ogy. In particular, the anti-screening predictions 
of f p4"p5| ] can be excluded. In addition, the exis- 
tence of a screening effect can be proved by he- 
lioseismology. 

REFERENCES 

1. S. Degl'Innocenti, W.A. Dzicmbowski, G. 

Fiorentini and B. Ricci, Astrop. Phys. 7 
(1997) 77. 



2. LOWL1+BISON: S. Basu et al., MNRAS 292 
(1997) 234; GOLF: S. Turck-Chieze et al, So- 
lar Phys. 175 (1997) 247; GONG: J.Harvey 
et al., Science 272 (1996) 1284; MDI: E.J. 
Rhodes et al. Solar Phys. 175 (1997) 287; 
LOWL1: S. Basu et al. Ap.J. 460 (1996) 
1064; GOLF2: Thierry et al. A&A 355 (2000) 
743. 

3. J.N. Bahcall, M.H. Pinsonneault, S. Basu and 
J. Christensen-Dalsgaard, Phys. Rev. Lett. 
78 (1997) 171. 

4. J.N. Bahcall, M.H . Pinsonneault and S. Basu, 
astro- ph/0010346! . 

5. A.S. Brun, S. Turck-Chieze and J. P. Zaha,Ap. 
J. 525 (1999) 1032. 

6. Fiorentini G. and Ricci B. Nucl. Phys. B. 
Proc. Suppl. 81 , 95 (2000). 

7. S. Degl'Innocenti, G. Fiorentini and M. Lis- 
sia, Nucl. Phys. B. Proc. Suppl. Series 43 
(1995) 66. 

8. V. Castellani, S. Degl'Innocenti, G. Fioren- 
tini, M. Lissia and B. Ricci, Phy. Rep. 281 
(1997) 309. 



7 



9. Borexino collaboration, C. Arpesella et al., 
"Borexino at Gran Sasso: proposal for a real 
time detector for low energy solar neutrinos" , 
internal report INFN, Milano, 1992. 

10. R.S. Raghavan, Phys. Rev. Lett. 78 (1997) 
3618. 

11. B. Ricci, F.L. Villante and M. Lissia, Phys. 
Lett. B 456 (1999) 214. 

12. GALLEX collaboration, XIX International 
Conference on Neutrino Physics and Astro- 
physics (Neutrino 2000), Sudbury, June 2000 
to appear on Nucl. Phys. B (Proc. Suppl.). 

13. SAGE collaboration, XIX International Con- 
ference on Neutrino Physics and Astrophysics 
(Neutrino 2000), Sudbury, June 2000 to ap- 
pear on Nucl. Phys. B (Proc. Suppl.). 

14. B. Cleveland et al., Ap. J. 496 (1998) 505. 

15. B. Ricci and F.L. Villante, Phys. Lett. B 488 
(2000) 123. 

16. V. Castellani, S. DegPInnocenti and G. 
Fiorentini, Astron. Astroph. 271 (1993) 601. 

17. J.N. Bahcall, Neutrino Astrophysics, Cam- 
bridge University Press, 1989 

18. Adelberget et al. Rev. Mod. Phys. 590 (1998) 
1265. see also : NACRE Collaboration, Nucl. 
Phys. AQ5Q (1999) 3. 

19. E.E. Salpeter, Australian J. Phys. 7, 373 
(1954). 

20. J.N. Bahcall et al. |astro-ph/0010055| (2000). 

21. G. Fiorentini, B. Ricci and F.L. Villante, 



astro-ph/0011130 



22. H.E. Mitler, Ap.J. 212, 513 (1977). 

23. G. Shaviv, Nucl. Phys. B Proc. Suppl, 38, 81 
(1995) 

24. V.N. Tsytovich, Astron. Astrophys. 356, L57 
(2000). 

25. V.N. Tsytovich and M. Bornatici, to appear 
on Plasma Phys. Rep. (2000). 



